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ABSTRACT 

We have computed a grid of high resolution synthetic spectra for cool stars 
(2500 < T eff < 6000 K) in the wavelength range AA 6000 - 10200 A, by 
employing an updated line list of atomic and molecular lines, together with 
state-of-the-art model atmospheres. 

As a by-product, by fitting TiO bandheads in spectra of well-known M 
giants, we have derived the electronic oscillator strengths of the TiO 7', 5, e 
and 4> systems. The derived oscillator strenghts for the 7', e and systems 
differ from the laboratory and ab initio values found in the literature, but are 
consistent with the model atmospheres and line lists employed, resulting in a 
good match to the observed spectra of M giants of known parameters. 

The behavior of TiO bands as a function of the stellar parameters T e ff, log g 
and [Fe/H] is presented and the use of TiO spectral indices in stellar population 
studies is discussed. 

Subject headings: Atomic and molecular data, Stars: atmospheres, M giants, 
globular clusters 
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1. INTRODUCTION 

A series of TiO spectral features are present in the integrated spectra of normal 
galaxies. Measurements of TiOi and Ti0 2 indices as defined in Burstein et al. (1984) and 
Worthey et al. (1994), and several indices given in Bica & Alloin (1986) are available in the 
literature. However, the use of TiO bands as stellar population indicators is not widespread. 

One clear characteristics of TiO bands is the dependence of their intensity on effective 
temperature. The effective temperature (T eff ) scale of M giants has been the subject of 
many studies in the past, both from the theoretical and observational standpoints - see the 
recent work by Bessell et al. (1998, hereafter BCP98) and Lejeune et al. (1998) The main 
source of uncertainty in the temperature scale of M giants ir related to uncertainties in the 
stellar mass, surface gravity and metallicity (Alvarez & Plez 1998). 

In this work, we use a revised atomic and molecular line list and state-of-the-art model 
atmospheres to synthesize the spectra of M giants in the near-infrared (NIR). A comparison 
of synthetic spectra with observations of cool red giants permits a calibration of T cff s for 
cool giants of the bulge metal-rich globular cluster NGC 6553. 

A grid of synthetic spectra in the wavelength range AA 6000 - 10200 A is computed 
and equivalent widths of TiO bands as a function of stellar parameters (T eff , log g, [Fe/H]) 
are measured. 

In Sec. 2 the observations are reported. In Sec. 3 the calculation of synthetic spectra 
is described. In Sec. 4 we derive effective temperatures for M giants of the globular cluster 
NGC 6553, based on the intensity of NIR TiO bands, and we derive electronic oscillator 
strengths for the TiO 7', S, e and systems. In Section 5 the computation of a grid of 
synthetic spectra in the NIR and measurements of equivalent widths of TiO features are 
presented; using the grid of synthetic spectra, the composite spectrum for a single-age 
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stellar population is built and compared to the observed integrated spectrum of the globular 
cluster NGC 6553. In Section 6 conclusions are drawn. 

2. OBSERVATIONS 

We have observed 6 stars among the coolest giants of the bulge metal-rich globular 
cluster NGC 6553 and two well-known field M giants HR 625 and HR 3816. The selection 
and identification of the stars in NGC 6553 were based on the images and Colour-Magnitude 
Diagrams (CMDs) given in Ortolani et al. (1990) and Guarnieri et al. (1998). The list of 
stars is shown in Table 1. 

The spectra of individual stars from the RGB tip of NGC 6553 were collected at the 
Cassegrain focus of the 1.6m telescope at the Laboratorio Nacional de Astrofisica (LNA), 
Pico dos Dias, Brazil, in 1997 June, using a EEV CCD of 770x1152 pixel (LNA CCD # 
48) of pixel size 22.5/ix22.5/im, and a grating of 8301/mm, with blaze angle at 7490 A. This 
combination gives a spectral resolution of AA ~ 3 A, centered at A8000 A with a spectral 
coverage of AA ~ 1350 A. In order to enlarge the spectral interval to include bands from all 
the relevant electronic transitions of TiO, we collected two spectra for each star, centered 
at A7570 and A8235 A with an overlap in the interval AA7570-8250 A. At the end of the 
reduction process, the spectra were combined and the final spectral interval is AA 6865-8905 
A. The flux calibration shows a systematic turnover around A 8700 A, which affects the 
spectrum shape at the < 5% level. The effect of this photometric inaccuracy on relative 
band intensities is negligible. 

An integrated spectrum of NGC 6553 and the M giant HR 625 were observed at LNA 
in 1996 July and in 1996 August respectively, at the Cassegrain focus of the 1.6m telescope, 
using a grating of 3001/mm, with blaze angle at A10000 A. Employing a 1024x1024 pixels 
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thin back-illuminated SITe CCD (LNA CCD # 101), of pixel size 24/zmx24/im, a spectral 
coverage AA 6020 - 10380 A, with an effective resolution of 8A, was obtained. 

The M giant HR 3816 was observed in 1996 February at the 1.5m telescope of the 
European Southern Observatory - ESO using the Ford Aerospace ESO CCD # 24, of 
2048x2048 pixels, with pixel size 15/xmxl5//m, together with a grating of 1200 1/mm (ESO 
grating # 11), giving a dispersion of 66 A/mm or a spectral resolution of AA ~ 2 A. 

Reductions were carried out using IRAF, following the usual steps: bias subtraction, 
flat-field and illumination correction, spectrum extraction and wavelength calibration. The 
spectra were flux calibrated through observations of three spectrophotometric standards per 
night, taken from Hamuy et al. (1994). Telluric lines have been removed through division 
by spectra of high vsini B stars. 



3. SPECTRUM SYNTHESIS 

A grid of synthetic spectra in the wavelength range AA 6000 - 10200 A is computed, 
using a revised version (Barbuy, B., Perrin, M.-N., in preparation) of the code described 
in Barbuy (1982), which incorporates the computation of molecular lines to the code for 
atomic lines by Spite (1967), assuming local thermodynamic equilibrium. The elemental 
abundances adopted are from Grevesse et al. (1996). The atomic lines are identified 
in the solar spectrum and oscillator strengths are obtained through a fit of synthetic 
spectra to the solar spectrum (Delbouille et al. 1973). The molecular systems taken 
into account are the CN (A 2 IT-X 2 £) red system, C 2 (A 3 r± 9 -X 3 r± u ) Swan system, TiO 
7 (A 3 $-X 3 A), 7 '(B 3 n-X 3 A), S^U-a^A), e(E 3 ri-X 3 A) and ^(b^-d 1 ^) systems, and the 
FeH (A 4 A-X 4 A) electronic system, which contains the Wing-Ford band. 

The detailed study of the C 2 Swan, CN red and TiO 7 systems and respective molecular 
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constants employed can be found in Barbuy (1985), Erdelyi-Mendes & Barbuy (1989), 
Barbuy et al. (1991), and Milone & Barbuy (1994). The list of FeH lines is described in 
Schiavon et al. (1997). The line lists and corresponding molecular constants for the TiO 7', 
5, e and systems, adopted from Jorgensen (1994), are discussed in Section 4. 

The photospheric models employed are from Kurucz (1992), Plez et al. (1992) and 
further unpublished models by Plez (1997), and Allard & Hauschildt (1995). The values 
of effective temperatures, surface gravities and metallicities adopted in the grid of models 
are reported in Table 4 (Sec. 5). Microturbulent velocities adopted were v t = 1.0km/ s 
and 2.5km/s for logg > 3.0 and logg < 2.5 respectively. The Plez et al. (1992) grid, 
which was employed in the computation of spectra in the range 2500 <T eff < 3400 K is 
only available for [Fe/H]=0.0, so that for [Fe/H] = -0.3 and -0.6 we have used the models 
for [Fe/H] =0.0, and the nonsolar metallicities were taken into account only in the input 
elemental abundances. 

The calculations were carried out in steps of 0.02A and the resulting spectra were 
rebinned to a step of lA and convolved with a gaussian function simulating an instrumental 
profile of FWHM=2A. We note that for the M giants studied here, it is important to take 
into account, in the calculation of continuum opacity, the H2 formation, which reduces the 
opacity due to HI and HII species. 

The photospheric C and N values are modified in giants due to convective mixing so 
that we adopted [C/Fe] = -0.2 and [N/Fe] = +0.4 for giants of log g < 2.0, and solar ratios 
for higher gravity stars. The change in C and N affects the TiO and CN bands: CN bands 
become stronger and, because of C deficiency, TiO bands also become stronger, since there 
is less CO association. 
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4. TiO BANDS AND THE T e// SCALE OF M GIANTS 

TiO bands are the main opacity source in the atmospheres of stars cooler than 4000K 
in the spectral region under study, and M giants are the dominant contributors to the 
integrated light of metal-rich and old composite systems such as bulge globular clusters and 
normal galaxies in the NIR. 

In Figure 1 are shown the contributions of the most important electronic systems of 
TiO in the studied spectral region, for a star of T e fr=3200K, logg=0.0 and [Fe/H]=-0.3. 
For A & 8000 A, the opacity is dominated by the 7 and 7' systems, whereas for A ^ 8000 
A, the 5 and e systems are dominant, the opacity due to the system being important only 
for A > 10000 A. 

The synthesis of the 7 system was already presented in Milone & Barbuy (1994): 
the laboratory line list by J.G. Phillips (see Phillips 1973) is used, and Honl-London and 
Franck-Condon factors were computed following the equations given by Kovacs (1979) and 
Bell et al. (1979) respectively; the electronic oscillator strength f d adopted here is given in 
Table 3 and discussed in Sec. 4.2. We have chosen to use the laboratory line list because it 
reproduces better the shape of the (0,0) 7 system bandhead (see discussion in Sec. 4.2.1). 

The line lists used for the computation of the 7', 5, e and (p systems, including 
the normalized product of the Franck-Condon and Honl-London factors, were the ones 
computed by Jorgensen (1994), which were kindly made available to us by R. Bell. For 
these systems, while the Honl-London and Franck-Condon factors were adopted directly 
from Jorgensen (1994), we derived eletronic oscillator strengths f e i values by fitting the 
main vibrational bands in the spectra of M giants. 

The fit to the intensity of TiO bandheads was applied to M stars of known stellar 
parameters (T eff , log g, [Fe/H], v t ): two well-known field M giants and M giants along the 
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cool red giant branch of the metal-rich globular cluster NGC 6553. We searched for VO 
bands in our spectra by following identifications by Barnbaum et al. (1996), and found 
no evidence for the presence of these bands, even for the coolest stars of our sample. The 
determination of the oscillator strengths proceeded in two steps, described in detail in the 
subsections below: (a) fine-tuning of the effective temperature of cool giants of NGC 6553, 
based on their TiO 7 system bandheads; (b) derivation of electronic oscillator strengths 
of the 7', e, 5 and systems, by fitting their main bandheads, adopting the temperatures 
established in (a). 



4.1. Stellar parameters for M giants of NGC 6553 

In a detailed analysis of individual stars of NGC 6553, Barbuy et al. (1997) obtained 
a metallicity [Fe/H] = -0.35 and a titanium-to-iron abundance ratio of [Ti/Fe] ~ +0.3. 
In a further analysis of these data, Barbuy et al. (1998) have confirmed that, in spite of 
a temperature uncertainty of 200 K, [Ti/H] « is found in all cases. As also discussed 
in Bruzual et al. (1997), a metallicity of [Fe/H] = -0.3 to -0.4 combined to enhanced 
[a-elements/Fe] = +0.3 to +0.4 result in Z/Z « 1 or [M/H] & 0. Furthermore, given 
that in our calculations the TiO bands are the dominant feature, and that [Ti/H] R3 is 
found for NGC 6553, it is clear that an overall metallicity [M/H] = 0.0 is to be used. In 
conclusion, we adopt [Fe/H] = 0.0 for NGC 6553, where by [Fe/H] we mean in fact [M/H]. 

For our sample of M giants, temperatures are first derived from the (V-K) colours 
given in Guarnieri et al. (1998). In Table 2 are given the observed (V-K), for which 
reddening corrections were applied using E(B-V) = 0.7 and E(V-K)/E(B-V) = 2.744 
(Rieke & Lebofsky 1985), and temperatures derived through the calibration of Ridgway et 
al. (1980), LCB98 and BCP98. 
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The 7 system bandheads were used to establish our T eff scale for the cool stars, given 
that among the TiO systems considered in this work, it presents the best settled molecular 
constants available. 

4-1.1. Fits to reference stars 

The fit to the field M giant HR 625, shown in Fig. 2, resulted in stellar parameters 
(Tea, log g, [Fe/H]) = (3600, 1.0, 0.0). This T cff is ~150K lower than the one estimated 
by Schiavon et al. (1997), on the basis of the spectral type M2 (Hoffleit & Warren, 1991) 
and the calibration by Fluks et al. (1994). In fact, the fit of the FeH bands in this star, 
adopting T eff = 3750 K was not very satisfactory in Schiavon et al. (1997). 

The fit to the vibrational sequence Av = (v',v") = [(0,0), (1,1), (2,2), (3,3), etc.] 
'composite bandhead' of the 7 system at A 7150 A to the spectra of stars in NGC 6553 
(Table 2), was applied in a fine-tuning of their effective temperatures, which are given in 
the column 8 of Table 2. For each T e ff we adopted a surface gravity based on the T e ff x log g 
relation given by LCB98, which is based on the evolutionary tracks of Schaller et al. (1992) 
for 1 M (column 9 of Table 2). 

The T e ff estimations given in Table 2 indicate that (i) temperatures from the (V-K) 
empirical fits of BCP98 are lower than those derived from the theoretical values of their 
Table 5; (ii) our TiO-based effective temperature estimations show a good agreement with 
the resulting temperatures from (V-K) calibrations by Ridgway et al. (1980), LCB98 and 
BCP98 (empirical). The T(TiO) values are adopted here. Taking into consideration the 
uncertainties in metallicity and log g, we estimate an error of AT c ff=150K for the stars 
hotter than T e s ~3200K, for which the gravity effect on band intensities is minimmum 
(see discussion on Sect. 5). For cooler stars, gravity effects are more important, so that 
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uncertainties can reach AT cfr =250K. Let us note that for T cff < 3200 K, the TiO 7 Av=0 
bandhead becomes saturated, so that for the star 291 (V12) the derived temperature is less 
precise - see also discussion on saturation of TiO bands in Sec. 5. 

In Fig. 3 the fit of a synthetic spectrum of the Av=0 vibrational sequence of the 7 
system to the star 288 of NGC 6553 (T cff = 3400 K) is shown. 

4.2. Electronic oscillator strengths for the TiO 7', e, 5 and systems 

The final f e i values obtained are given in Table 3, together with other determinations 
collected from the literature: the values by Davis et al. (1986) were obtained from 
correlations between equivalent widths of rotational lines and lifetime estimates for the a 
and (3 systems; Brett (1990) obtained semi-empirical values through the fit of observed 
spectra of M giants; Jorgensen (1994) reports f ei values used in his calculations; Hedgecock 
et al. (1995) derived / e ; from lifetime measurements; Langhoff (1997) ones were obtained 
from ab initio calculations; Alvarez & Plez (1998) present updated f e i values from recent 
laboratory measurements or calculations, which they also checked through comparisons 
between observed and calculated spectra of M giants. 

The f e i values used in our calculations were derived in this work as follows: 

(a) Adopting the T e fjs as described in Sec. 4.1, we derive f e i values for the e and 5 
systems, by fitting the observed intensity of their Av=0 vibrational bands, as shown in Fig. 
4. In the case of the e system, we had to shift the wavelengths of all lines by AA = +76.0 
A in order to correctly reproduce the position of the (0,0) bandhead at « A8400 A. This 
Ao shift is compatible with an uncertainty of the order of 100 cm^ 1 in the electronic term 
T c of the upper level (E 3 II) of the transition (see Jorgensen 1994; Huber & Herzberg 1979; 
Simard & Hackett 1991). We find T c = 11829.3 for the E 3 II level, by shifting Jorgensen's 
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vable. The value obtained for this system (f eJ = 0.07) is in disagreement with all values 
in the literature by a factor > 10. Recent laboratory determinations by Lundevall (1998) 
give a lifetime of 4.0±0.2 /j,s for the E 3 I1 level, wherefrom a f e ; even lower than Langhoff's 
would be derived. Therefore we preferred to omit our result for the e system in Table 3. 
For our purposes of computing a library of synthetic spectra, it is only important that the 
combination of line list, Franck-Condon factors, Honl-London factors and i e \ match the 
observed spectra. We suggest that the line list of Jorgensen (1994) for the e system is to be 
revised. 

(b) The intensity of the vibrational sequences given by Av = [(v',v") = (0,0), (1,1), 
(2,2), (3,3), etc.] and Av = -1 [(v',v") = (0,l), (1,2), (2,3), etc.] of the i system, at « 
A6200 and A6700 A respectively, were used to obtain its electronic oscillator strength f eh 
by fitting the spectrum of HR 625. 

(c) The f e i of the system was adopted from the ratio of the theoretical f e is derived 
by Langhoff (1997) for the 5 and systems, giving ffj ' f 8 d ~ 0.2 . 

We adopted the procedure of fitting the bands from the different systems to the 
spectrum of the star NGC6553-288, for which the T e g determined from the intensity of 
the bands from the 7 system is in best agreement with the T e ffS derived from the various 
colours. The f e ;s so derived for the other systems resulted in good fits to the spectra of all 
the remaining stars, of Table 2, adopting T(TiO). 

The values derived in this work are in most cases similar to the average of the values 
estimated by other authors. The agreement between different authors is best for the 5 
system. For the 7 system, our adopted value of f e / = 0.12 corresponds essentially to a mean 
of the laboratory values U = 0.165 and 0.092 by Davis et al. (1986) and Langhoff (1997). 



- 12 - 



4-2.1. Errors in the electronic oscillator strenghts 

The main sources of uncertainty in the derivation of our astrophysical f e ;s are the 
uncertainties in T cff , logg, Ti and O abundances and the adopted microturbulent velocity 
(v t ). Another potential source of uncertainties is a possible incompleteness in the laboratory 
line list adopted for the 7 system. The latter is an important issue given that the f ei s of all 
the remaining systems are tied to the T e gfS derived from the intensities of the bands of the 
7 system. 

We have adopted a fixed v t = 2.5km/ s for giants. For v t = 1.5km/ s there would occur 
an increase in the estimated f e i of 15%. 

The uncertainty in the C and N abundances is estimated by computing the 7 band 
with [C/Fe] = [N/Fe] = 0.0: the change relative to the adopted values of [C/Fe] = -0.2 and 
[N/Fe] = 0.4 is of the order of 2% in the estimated f e i value. 

Adopting either a T e g- or a metallicity differing in 200K and 0.3dex respectively 
results in a shift of 20% in the resulting f e ;s. As the log g sensitivity of the TiO bands at 
T cff ~3400K is very low (see Sec. 5), the influence of this parameter upon the derived f e ;S is 
very weak. An error of 0.5dex in log g results in a difference of only 4% in the resulting 
f e ;S. Taking into consideration all the above sources, the final uncertainty in our i e i values 
is around 30%. 

We also computed synthetic spectra using the theoretical line list from Jorgensen 
(1994) for the 7 system, in order to check for a possible incompleteness of the Phillips 
laboratory line list. As a result, the f e « that should be adopted with the Jorgensen line list 
is only marginally lower (8%) than the one that is compatible with the Phillips line list. 
Taking into consideration the other sources of uncertainty discussed above, we conclude 
that such an incompleteness in the Phillips line list has only a mild effect in our analysis. 
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Moreover, we chose to adopt the Phillips line list because it gives a better reproduction of 
the band shapes. The Jorgensen line list gives too high opacities at high J and too low 
opacities near the bandheads, so that it is very difficult to correctly reproduce the whole 
observed band shapes when employing this line list. 

The use of the f e ; value by Langhoff (1997) or Hedgecock et al. (1995), which are lower 
than ours by respectively 30% and 50%, would imply too low temperatures in comparison 
to those derived from the photometry (see below). We consider that these values are 
incompatible with our data. 

5. A GRID OF SYNTHETIC SPECTRA: TiO vs. STELLAR 

PARAMETERS 

We have computed a grid of synthetic spectra for the stellar parameters indicated in 
Table 4. A list of indices which measure TiO bands, as detailed in Table 5, was measured 
for the whole grid of spectra, in order to study their behavior as a function of stellar 
parameters. 

In Fig. 5 are shown the plots of equivalent widths of TiO bands measured (cf. Table 
5), for different log g and for a fixed metallicity [Fe/H] = -0.3, as a function of T e fj. The 
bandheads at A 6000, 6600, 7150, 7600 and 8400 A show similar behaviors, as expected. 
The dependence on gravity is also shown. In Fig. 6 the behavior of the indices, for a fixed 
logg = 0.0, and different [Fe/H], as a function of T c fr, shows the strong dependence on 
metallicity. 

The Ti02 index as defined in Burstein et al. (1984) (Figs. 5b and 6b) is much less 
sensitive to surface gravity than EWeooo- In order to explain this, we plot in Fig. 7 the 
synthetic spectra for Teff = 3400 K, log g = 0.5 and [Fe/H] = 0.0, where the Ti0 2 index as 
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defined by Worthey et al. (1994) (the index is slightly modified relative to Burstein et al. 
1984), as well as our EW 60 oo index are indicated (see also Table 5). It can be seen that the 
Ti02 index measures the bottom of the TiO Av=0 vibrational sequence of the 7' system, 
where line saturation is stronger. The index is clearly saturated for T eff < 3500 K (for 
[Fe/H] = 0.0), so that its sensitivity to both T eff and \ogg is very small in this T eff range - 
see Fig. 6b. 

Figs. 5a-f and 6a-f show that the redder the TiO feature, two effects occur: the 
threshold of detection of the bands starts at progressively lower temperatures, and the 
bandheads saturate for cooler temperatures. For example, EW 60 oo and EW 8 4oo rise m 
strength for T eff > 4000 and 3400 K whereas saturation occurs at T cff < 3000 and < 2750 
K respectively. It has to be noted that in all plots of Figs. 5 and 6 the equivalent width 
corresponding to T eff =2500 K is underestimated, due to the use of log g = -0.5, instead of 
log g = -1.0, more appropriate for this temperature, since no models are available for log g 
= -1.0 in the Plez grid. 

One conclusion that can be drawn from these plots is that for composite stellar 
populations of different mean temperatures, the use of different indices among those listed 
in Table 5 will be more appropriate in each case. For a very metal-rich and cooler stellar 
population of a massive elliptical galaxy, EW8400 would be more sensitive to changes in line 
strength, whereas for a metal-poor globular cluster the bluer indices are more suitable. 

Fig. 8 shows the Av=+1 vibrational sequence of the S system, fitted to HR 3816. The 
features are clearly visible, although much weaker than the bluer ones, even for such a low 
temperature as T cff = 2750 K. The behavior of the corresponding index EW 98 oo vs. T cff 
(Fig. 9) indicates that this band essentially does not saturate in the temperature range 
considered. 
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5.1. Integrated spectrum of NGC 6553 

We concentrate our attention on the bulge cluster NGC 6553, for which we have a 
wealth of information: i) V and I photometry from the Hubble Space Telescope (Ortolani 
et al. 1995); it) J and K photometry from the ESO 2.2m telescope (Guarnieri et al. 1998), 
resulting in reliable determinations of reddening and distance to the cluster; in) detailed 
abundance analyses of stellar members (Barbuy et al. 1998), giving [Fe/H]~-0.4 and 
[a/Fe]~+0.4; iv) medium-resolution spectra of stars from the tip of the red giant branch 
(Sec. 2) in the NIR, including several TiO bands, and v) the integrated cluster spectrum 
covering; the spectral interval AA 6000 - 9700 A. 

From the observed CMD, it is possible to compute the integrated synthetic spectrum 
of NGC 6553 in a straightforward manner. Such method was employed by de Souza et al. 
(1993), Barbuy (1994) and Milone et al. (1995) where synthetic spectra were used, and 
Santos et al. (1995) where the stellar library of Jacoby et al. (1984) was adopted. On the 
other hand, the integrated spectrum of NGC 6553 was reproduced in Bruzual et al. (1997), 
by using isochrones and several spectral libraries (synthetic and observed). 

The CMD of NGC 6553 was transformed into a loggxT eff diagram, using the 
T ef f x loggx[Fe/H]x(V-I) relations from BCP98 for all the stars, except the ones at the 
RGB tip. For the latter, we adopted the T e ff scale derived in Section 4. For the stars 
below the turn-off, we adopted a Salpeter IMF plus the M-dwarf theoretical main sequence 
from Baraffe et al. (1995). Each star in the loggxT eff plane has its absolute I magnitude 
(computed from distance and reddening estimated by Guarnieri et al. 1998), and a 
corresponding synthetic spectrum, which is chosen to be the one with closest atmospheric 
parameters. The integrated synthetic spectrum is then obtained by summing up all the 
selected individual stellar spectra, weighted by M/. In Fig. 10 the result is compared with 
the observed integrated spectrum of NGC 6553, adopting [M/H] = 0.0. 
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From the above calculations, we can estimate the contribution of stars from all 
evolutionary stages to the integrated light of the cluster. This is shown in Fig. 11, where 
it can be seen that the few stars in the RGB tip largely dominate the cluster light in the 
spectral region under consideration. 

We note that the implementation of this grid of synthetic spectra to the code of 
evolutionary populations synthesis by Bruzual & Chariot (1998) is under way. 

6. CONCLUSIONS 

We present a new grid of synthetic spectra in the interval AA 6000-10200A. This grid 
is computed by employing state of the art model atmospheres and has a higher resolution 
than the ones currently in use in stellar population synthesis in the same spectral interval. 

A fit to NIR spectra of stars from the RGB tip of NGC 6553 enabled us to derive a 
T e ff scale for the M giants, based on the TiO 7 system band intensities, on the one hand, 
and also to determine a set of semi-empirical electronic oscillator strengths for the 7', 5, e 
and 4> systems of the TiO molecule. 

The adequacy of our spectral library for stellar population synthesis purposes was 
tested for the case of NGC 6553, a well studied bulge globular cluster, for which detailed 
abundance analyses of individual stars and well-defined CMDs are available. In this 
framework, we obtained a good match to the ISED of this cluster. This result warrants the 
use of our spectral grid in stellar population synthesis studies in the NIR. 

The authors are indebted to Robert Kurucz, France Allard and Bertrand Plez for 
providing their model atmospheres, Roger Bell for the theoretical TiO line list and for 
helpful comments, and to Michael Bessell for making available his results in advance of 
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Table 1. Log of the observations of field and cluster M giants. 



star 


V 


«1950 


"1950 


date 


exp.(m) 


AA (A ) 


HR 625 


6.10 


02 h 06 m 23.4 s 


-18°00'55.5" 


05.08.96 


2 


5900-9200 


HR 3816 


6.10 


09 30 59.2 


-62 34 01 


01.02.96 


0.5 


8090-10040 


NGC 6553-int. 




18 05 11 


-25 55 06 


19.07.96 


30 


5900-9200 


NGC 6553-291 


16.92 






10.06.97 


30 


6865-8251 


NGC 6553-125 


15.89 










6865-8905 


NGC 6553-288 


15.91 










ii 


NGC 6553-47 


15.35 












NGC 6553-118 


15.70 










ii 


NGC 6553-211 


15.61 










ii 
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Table 2. Frame coordinates, observed (V-K) colours and temperatures of 6 cool giants of 
NGC 6553. Coordinates correspond to the frame in K band shown in Fig. 1 by Guarnieri 
et al. (1998); star numbers and colours are also from Guarnieri et al. The field stars HR 
625 and HR 3816 are also included in the Table. Relations colour-temperature by Bessell et 
al. (1998) - BCP (theoretical/empirical), Lejeune et al. (1998) - LCB98, and Ridgway et 
al. (1980) are used. Gravities adopted are given in the last column. 



star 


X 


Y 


(V-K) 






T 

± Ridgway 




log g 


291/V12 


-4.7 


324.0 


9.75 


3250/3075 


3100 


3060 


3200 


0.0 


125 


30.1 


276.5 


7.21 


3520/3476 


3550 


3530 


3500 


0.7 


288 


10.8 


334.0 


7.93 


3420/3367 


3400 


3410 


3400 


0.5 


47 


63.7 


286.2 


6.96 


3550/3518 


3590 


3570 


3400 


0.5 


118 


30.0 


281.0 


7.15 


3530/3476 


3550 


3540 


3500 


0.7 


211 


14.1 


224.9 


3.89* 








3500 


0.7 


HR 625 






4.48 


3660/3670 


3700 


3680 


3600 


1.0 


HR 3816 














2750 


-0.5 



Note. - * For the star 211 the (V-K) colour appears to be in error, since the 
spectrum is incompatible with a high temperature 
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Table 3. Electronic oscillator strengths from the literature (Davis, Littleton & Phillips 
1986 - DLP86; Brett 1990 - B90; Jorgensen 1994 - J94; Hedgecock, Naulin & Costes 1995 - 
HNC95; Langhoff 1997 - L97; Alvarez & Plez 1998 - AP98) compared to the value found in 

the present work 



system 


DLP86 


B90 


J94 


HNC95 


L97 


AP98 


This Work 


7' 


0.138 




0.14 


0.093 


0.108 


0.0935 


0.06 


7 


0.165 


0.22 


0.15 


0.079 


0.092 


0.0786 


0.12 


5 


0.05 


0.05 


0.048 




0.096 


0.048 


0.05 


e 




0.006 


0.014 


<0.0056 


0.002 


0.0023 







0.05 


0.05 


0.052 




0.018 


0.0178 


0.01 



Table 4. Stellar parameters and corresponding photospheric models used to build the grid 
of synthetic spectra; models are from Kurucz (1992) - K92, Plez (1997) - P97, Plez et al. 
(1992) - PBN92 and Allard & Hauschildt (1995) - AH95 



T cff /step log g/step [Fe/H] models 



4500- 


-6000/250 


0.0-5.0/0.5 


-0.3,0.0 


K92 


4000- 


-4500/250 


2.0-5.0/0.5 
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PBN92 


2750- 


-3000/250 


-0.5-1.5/0.5 


-0.6,-0.3,0.0 




2500 




-0.5 


-0.6,-0.3,0.0 


» 


2800- 


-3200/200 


3.5-5.0/0.5 


-0.5,0.0 


AH95 


3300- 


-3700/200 


3.5-5.0/0.5 


-0.5,0.0 




2500 




5.5 


-0.5,0.0 
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Table 5. Definition of spectral indices measured 



Index 


T e// range 


Blue continuum 


Bandpass 


Red continuum 


EW6000 


< 


3600 


6143.95-6146.41 


6145.18-6540.93 


6539.58-6542.27 


EW6000 


> 


3600 


6067.47-6074.81 


idem 


idem 


J^VV 6600 


< 


3600 


6oo9.oo-6o4z.z ( 


p. c a n no 7n/io on 
654U.93- / U4o.2z 


(o4e. fO- food. 14 


EW6600 


> 


3600 


idem 


idem 


7594.40-7598.00 


EW7100 


< 


3600 


7046.36-7050.08 


7048.22-7550.92 


7548.70-7553.14 


EW7100 


> 


3600 


6539.58-6542.27 


idem 


7594.40-7598.00 


EW7 6 oo 


< 


3600 


7548.70-7553.14 


7550.92-8190.40 


8188.86-8191.93 


EW7 6 oo 
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3600 


7593.64-7598.33 


idem 


8259.46-8261.50 


EWs400 


< 


3600 


8188.86-8191.93 


8190.40-8853.40 


8850.67-8856.13 


EW8400 


> 


3600 


8259.46-8261.50 


idem 


idem 


EW9800 






9721.67-9725.00 


9726.00-9972.00 


9975.23-9981.52 


Ti02(Burstein) 






6069.00-6142.75 


6192.00-6273.25 


6375.00-6416.25 


Ti0 2 (Worthey) 






6068.375-6143.375 


6191.375-6273.375 


6374.375-6416.875 
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FIGURE CAPTIONS 

Fig. I. — Strongest TiO bands in the near infrared, corresponding to the 7, 7', e and 5 
electronic systems. 

Fig. 2. — 7 and 7' bandheads for HR 625. Solid lines: synthetic spectrum computed for 
T eff =3600 K, log g = 1.0, [Fe/H] = 0.0; dashed line: observed spectrum 

Fig. 3. — Av=0 and Av=-1 vibrational sequences of the 7 system for the star NGC 6553:288. 
Solid line: synthetic spectrum computed with T ef j=3400 K, log g = 0.5, [Fe/H] = 0.0; dashed 
line: observed spectrum. 

Fig. 4. — Av=2 vibrational sequence of the 7 system, Av=0 of the e system and Av=0 of 
the S system for NGC 6553:288. Solid line: synthetic spectrum computed with T eff =3400 
K, log g = 0.5, [Fe/H] = 0.0; dashed line: observed spectrum. 

Fig. 5.— TiO features vs. T cff for [Fe/H] = -0.3 and log g = -0.5, 0.0, 0.5, 1.0 and 1.5. 

Fig. 6.— TiO features vs. T cff for log g = 1.5 and [Fe/H] = -0.6, -0.3 and 0.0. 

Fig. 7.— Synthetic spectra for Teff = 3400 K, log g = 0.5 and [Fe/H] = 0.0, where the Ti0 2 
index as defined according to Worthey et al. (1994), and our EW 60 oo index are indicated. 

Fig. 8. — Av = +1 vibrational sequence of the 5 system in HR 3816. Solid line: synthetic 
spectrum computed for T cff =2750 K, log g = -0.5, [Fe/H] = 0.0; dashed line: observed 
spectrum. 

Fig. 9. — TiO feature EW 980 o vs. T cff . Symbols are as in Fig. 5 and open stars correspond 
to dwarfs. 

Fig. 10. — Integrated spectrum of NGC 6553. Dashed line: observed spectrum; solid line: 
synthetic spectrum computed as described in Sec. 5.1. 
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Fig. 11. — Contribution of stars in different evolutionary stages to the integrated light of 
NGC 6553. 
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